One-dimensional optical lattices and impenetrable bosons 
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We study the limit of large onsite repulsion of the one-dimensional Bose-Hubbard model at low 
densities, and derive a strong-coupling effective Hamiltonian. By taking the lattice parameter to 
zero, the Hamiltonian becomes a continuum model of fermions with attractive interactions. The 
leading corrections to the internal energy of a hard-core-boson (Tonks) gas as well as the (finite 
temperature) pair correlations of a strongly interacting Bose gas are calculated. We explore the 
possibility of realizing, in an optical lattice, a Luttinger liquid with stronger density correlations 
than the Tonks gas. A quantum phase transition to a charge-density- wave Mott insulator is also 
discussed. 
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The field of ultra-cold atoms is attracting much atten- 
tion in recent times. This is partly due to the possibility 
(in some cases already experimentally demonstrated) of 
realizing exotic quantum phase transitions 0. One of 
the most interesting examples is the transition between 
a superfluid and a Mott insulator recently observed in an 
optical lattice 0] • The prototypical model for this transi- 
tion is the Bose-Hubbard model, which in one-dimension 
takes the form: 



H B H = ~- [ b tn+An + h.< 



- Y X! n m (n m - 1), 

(i) 

where [frm,^/] = Sm,m'i and commute otherwise; n m = 
b\ n b m counts the number of bosons at the m-th site, 
where m = 1, . . . , M a (M a is the number of lattice sites). 

As it has been recently discussed by Biichler et al. Q 
(see also HE)), in one-dimension the superfluid-Mott 
insulator transition can be described by the sine-Gordon 
model: 



H SG = H gauss [6, 0, K]+g / dx cos [28(x) + 2irx5] . (2) 
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where 



H = ^ / dx 



{d x ef 

K 



K{d x(l 



(3) 



describes the (long wave-length) gaussian fluctuations of 
the phase <j>{x) and density 9{x) fields & = P° ~ a ~ 1 
is the mismatch between the mean density, p Q , and the in- 
verse lattice parameter a -1 . The sine-Gordon model ex- 
hibits a Kosterlitz-Thouless (KT) transition [3,3 close to 
K = 2 for i5 = 0, and a commensurate-incommensurate 
transition for K < 2 as 5 is varied Q ■ The dimensionless 
parameter K is proportional to the compressibility |lOj| . 
and controls the spatial decay of phase and density corre- 
lations. Since in one dimension Bose-condensation can- 
not take place even at zero temperature, correlations de- 
cay at least algebraically. Thus, for a sufficiently large 
system at T = 0, the one-body density matrix displays a 



power- law behavior HQ: (*t( x )$(o)V ^ -M- 1 ^. Den- 
sity correlations, however, decay as 0, (p{ x )p(^)) ~ 
x~ 2 + A cos(2Trp )\x\~ 2K , where A is a constant. For a 
zero-range interaction, v(x) — g5(x), Lieb 0] found the 
exact eigenstates and eigenvalues of the Hamiltonian us- 
ing the Bethe ansatz. This allows to obtain K for every 
value of the dimensionless parameter 7 = Mg/h 2 p a (M 
being the boson mass). However, for this model K > 1 
always, and takes the smallest possible value (K = 1) as 
7 — > +00 and the system becomes a Tonks gas [jJ^J. In 
this limit, many of the long-distance correlations are the 
same as for a non- interacting Fermi system 0, 0, ^| . 

In this paper we will consider the large U limit of the 
Bose-Hubbard model, Eq. when the number of par- 
ticles, N , is less than the number of lattice sites. We 
show the following: 

1. For U 3> J, at low temperatures, Eq. can be 
effectively replaced by a model of spinless fermions 
with correlated hopping and attractive interac- 
tions: 
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where {c m ,c| n ,} = 5 m ,m', anti-commuting other- 



wise, and n 



2. When the lattice parameter a — > and the Bose- 
Hubbard model becomes a continuum model of 
bosons interacting with zero-range interactions (the 
so-called Lieb-Liniger model), Hp becomes |l9j : 



Hp = / dx 

2M J Q 

/ dxd x ^ F {x)d^ F {x)^ F {x)* F {x), (5) 
* Jo 

with the identifications typ(x m = ma) = ar x l 2 c m , 
L = M Q a, g F = 4h 4 /M 2 g, and fi/M = Jo 2 , 



2 



g = Ua. Using this effective fermionic Hamil- 
tonian we provide elementary derivations for the 
0(7 _1 ) correction to internal energy as well as for 
the (temperature-dependent) pair correlation func- 
tion: g 2 = Po 2 {^(x)^(x)^{x)^{x)), where ^(x) 
is the boson field operator. 

3. By including the interaction with the nearest neigh- 
bors, we argue that it is possible to realize in an 
optical lattice a Luttinger liquid with K < 1. In 
this regime a boson system resembles a system of 
fermions with repulsive interactions. In particular, 
the system can undergo a quantum phase transition 
to a charge-density-wave Mott insulator in which 
the boson density is modulated with twice the pe- 
riod of the optical lattice. 

We now proceed to derive the results listed above. Let 
us consider a one-dimensional optical lattice. As usual, 
this presumes the existence of a tight confinement in the 
transverse direction; the effect of the confinement in the 
longitudinal direction will be discussed at the end. As- 
suming that U is the largest energy scale in the prob- 
lem, i.e. U ~> max{J, T}, only states with site occu- 
pancy equal or less than one (i.e. n m < 1) will be avail- 
able because they have zero interaction energy. Thus we 
can split the Hilbert space into two complementary sub- 
spaces: H = Hp © Hq. We denote by P the projector 
onto the subspace Hp of states which obey the constraint 
ri m < 1 for all m; Q = 1 — P projects onto the comple- 
mentary subspace, Hq. Our goal in what follows is to 
find the effective Hamiltonian within the subspace Hp. 
To t his p urpose we use the perturbative expansion (see 
e.g. P): 



H cS = PH BH P - PH BH Q 



1 
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(6) 



Here we shall not need the terms denoted by "• ■ •" To 
0(J 2 /U) we obtain: 



~w E p ( b Ui + bl^)b m Qbi(b m+1 + b m ^)p. 

(7) 



m— 1 



The first term in the previous expression describes the 
hopping of impenetrable bosons. The projection opera- 
tors can be removed if we represent the boson operators 
by Pauli matrices: (fct, + b rn ) -> cr* , i(b m - b] n ) -> of„, 
and 2n m — 1 — > a^. Alternatively, as we shall do hence- 
forth, one can use Jordan- Wigner fermions (see e.g. pj) 
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b m = exp [inJ^Kni n l] c m, and n m = 6^& m 
ter restricting ourselves to Hp. In terms of the fermionic 
operators, 
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(8) 



We now turn our attention to the second term in 
Eq. |0. The hopping terms on the left side acting on 
P move particles from the neighboring sites to the m-th 
site. Therefore, in the resulting state, the m-th site can 
be either singly our doubly occupied. The former possi- 
bility is projected out by the operator Q. Thus we can 
replace Q by (n m — 1), which does nothing when n m — 2 
and vanishes when n m — 1. After some manipulations, 
we obtain: 
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rn — 1 



P(9) 

Finally, we perform the Jordan- Wigner transformation to 
remove the projector P |20|. Collecting all terms yields 
the Hamiltonian Hp given in Eq. This proves the 
first result announced above. 

Next we consider the continuum limit of Hbh and 
Hp. In both cases, if interactions are neglected, the 
Hamiltonian is diagonalized by Bloch waves: b p — 

M a 1 ^ 2 J2m=i e ~ tpXm bm (an identical expression applies 
to fermions, with c's replacing fc's). On the lattice, 
non-interacting particles disperse according to e(p) = 
— Jcos(pa). Taking a — * and expanding e(p) about 
p = allows to identify, by comparing with the free- 
particle dispersion e (p) = h 2 p 2 /2M, the coefficients of 
the quadratic term: Ja 2 = fr 2 /M. Therefore, by keep- 
ing the product J a 2 constant as a — > 0, one recov- 
ers the continuum model of free particles. In a simi- 
lar way, after introducing the continuum field operators 
^(x m ) — CL~ 1 ^ 2 b r rn the interaction term of Hbh becomes 
H int = § f Q L dx ^(x)^(x)^(x)^(x), where g = Ua and 
L = M Q a are kept constant as a — » 0. Let us now take 
the a — > limit of the interaction term in Hp. Using 
$>F(x m ± a) = a _1 / 2 c„ l± i k ^ F {x m ) ± ad x ^ F (x m ) + 



a 2 d 2 ^F(x m )/2 + 0(a J ), where x„ 



and intro- 



ducing this expression and its hermitian conjugate into 
Eq. (QJ, Eq. (J5J is obtained. To close this discussion, it 
is important to point out that the continuum limit which 
we have taken does not correspond to the usual contin- 
uum limit for fermions (see e.g. |1J) where one linearizes 
the dispersion around the Fermi points ±p F = ±7rp . 
Instead, we have expanded around the bottom of the 
band. This is in agreement with the way in which the 
Lieb-Liniger model is obtained from the Bose-Hubbard 
Hamiltonian pi) . 

To obtain the 0(7~ 1 ) correction to the internal en- 
ergy of a strongly interacting Bose gas, it is convenient 
to express Hp in Fourier components (henceforth we use 
periodic boundary conditions before taking the thermo- 
dynamic limit L — ► +oo): 

Hp=Y^ <P)4 C P + f£ E( P + ^ k C l+q C k-q C kCp. (10) 
p pkq 

The free energy is obtained from the logarithm of the 
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partition function (/? = 1/T): 



Z = Tr e 



-PiH-nN) 



Z (Texp 



ft/3 



da flint (cr) 



>o, 
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where is the Hamiltonian of the Lieb-Liniger model 
and AT the total particle number operator. In the second 
expression on the left hand-side Z Q = Tr p ((3,p,), and 



(•••)o means average with p Q (j3,p) 



-P(H a -u.N) 



/Z , 



i.e. over the non-interacting grand-canonical ensemble. 
The key to the calculations that follow is to realize that 
for 7 > 1, we can replace H by Hp. Thus, expanding 
Eq. in powers of gF to O(gp), 
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(12) 



To obtain the last expression we have used 
Wick's theorem, denoting the (time-ordered) 
free-fermion Matsubara-Green's function as 
G(x,a) = -£p<T[cp((7)4(0)]) o e<**/L; in the 
same expression the limit e — > _ must be taken. The 
internal energy can be obtained from the above result 
by using the standard expression U = — 9 gf Z + p, (N) , 
which in the present case yields: 



U 
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(13) 



where U — ^2 p c{p){n p ) is the internal energy of a 
free Fermi gas (which equals that of a Tonks gas), and 

(n p ) = (e^^)-^) + 1) _1 the Permi-Dirac distribu- 
tion. Using the Sommerfeld expansion |l2j, U /N = 
(h 2 p 2 /2M)(TT 2 /3)[l + T 2 /47r 2 + 0(T 4 )], where r = 
T/Td and T c i = h 2 p 2 /2M the degeneracy temperature in 
units where fee = 1- Introducing this result into Eq. (I13|) . 
we obtain for r 1 and 7 > 1: 



U(T) h 2 p 2 
N 2M 
4 
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(14) 



At zero temperature, this expression reduces to E gs = 
U(T = 0) = N (h 2 p 2 /2M)(n 2 /3)_(l - 4/7 + O^ 2 )), 
which agrees to leading order in 7 1 with the result ob- 
tained by Lieb and Liniger in Ref. 0. 

Equation i|12|l also allows us to obtain asymptotic ex- 
pressions for the pair-correlation function. As recently 
pointed out by Kheruntsyan et al. for a homoge- 

neous system this function can be obtained using the 
Hellmann-Feynman theorem: 



1 dlogZ 
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Applying this to Eq. l(T2*)l yields (e — > ) : 



.92 



l 2 Pi 



(5 x G(0,e)) 2 -G(0,e)92 G (0, e ) 



(16) 



This provides an explicit derivation of the expression 
employed in Ref. [13J, where Eq. I|16f) was used after ex- 
trapolating to finite temperatures the T = result ob- 
tained from the Bethe ansatz wave-function. We be- 
lieve that the present derivation is simpler, and does not 
require to deal with the Bethe ansatz solution. Using 
this expression we obtain in terms of 7 and r = T/T^ 
that g 2 (-f, t) = I (7r/7) 2 [l + (r 2 /47r 2 )] , in the quan- 
tum degenerate regime t«1, and .92(77"?") = 2r/ 7 2 for 
1 t <§C r ) 2 , where the gas is non-degenerate but bosons 
are still strongly interacting. 

Let us finally discuss some of the phenomena that char- 
acterize the regime where K < 1. We first discuss, very 
briefly (further details will be provided in a future pub- 
lication ^3), now to enter this regime in an optical lat- 
tice. In realistic situations, when U is the largest en- 
ergy scale of the problem, the description provided by 
Eq. is not complete. There is always present a small 
repulsion between bosons in neighboring sites 0, IbH ] 
which is described by the operator V n m n m +i, with 
< V <C U. When this is taken into account, the rele- 
vant Hamiltonian at large U becomes: 
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Jn 
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m ■[■ 



V ^2 n mn m +i + O 



U 2 



(17) 



where V' — V — J 2 /U. As we argue below, the K < 1 
regime can be reached by increasing V. To this purpose, 
we notice that both V and U are proportional to the 
strength of the interaction between the atoms, which can 
be tuned to larger values using Fcshbach resonances. For 
instance, in an optical lattice near half-filling (i.e. for 
N a w M a /2) and for sufficiently larg e U, the model of 
Eq. (fP7|l can be well approximated |lj| by the Heisenberg- 
Ising chain 15]. For this model, the expression for K fig 
obtained from its Bcthe-ansatz solution, 



K 



27r-2cos- 1 (U / /J)' 



(18) 



tells us that K < 1 can be achieved by making < 
V < J. One immediate consequence of the existence 
of a strong repulsion with nearest neighbors is the en- 
hancement of density fluctuations. For instance, if the 
system is tuned to a Luttinger liquid with K < 1/2, the 
density correlation function (which behaves as \x\ 



-2K 



at 



T = and L 



-00) will decay more slowly than the 



phase correlations (behaving as 



- 1 ' 2K ). This con- 



clusion remains true when the system is confined in a 
trap [lOj. Furthermore, in a trap the density profile ex- 
hibits Friedel oscillations of period equal to p" 1 near its 
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boundaries jToL . The amplitude of these oscillations 
is strongly suppressed for K > 1, and strongly enhanced 
for K < 1 . The higher visibility of the Friedel oscillations 
can provide an experimental signature that this regime 
has been reached. By contrast, in the current situation 
many of the experimentally available quasi-ld systems 
are quasi-condensates (K > 1) with a fluctuating phase 
but strongly suppressed density fluctuations. 

Furthermore, near half-filling N a ss M a /2 and for 
V > J, a quantum phase transition to a new type 
of (boson) Mott insulator can take place. The tran- 
sition is again of Kosterlitz-Thouless type and can be 
described by a sine-Gordon model of the form H = 



#gauss + g J dx cos [A9(x) + AttxS'] (see e.g. 0,0)- In the 
insulating ground state (V > J), the lattice occupancy 
n m follows the patterns 010101010 • • ■ or 101010101 ■ ■ •, 
i.e. the boson density is modulated with twice the lat- 
tice period. This could be observed by means of Bragg 
scattering. The excitations are solitons, interpolating be- 
tween the two degenerate ground states: For example 
• ■ • 010101101010 • ■ • This excitation is fractional and car- 
ries half a boson [TM Il8| . A longitudinal confinement can 
be treated as described in Ref. 
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